Scattering of charge carriers by point defects in bilayer graphene 
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Theory of scattering of massive chiral fermions in bilayer graphene by radial symmetric potential 
is developed. It is shown that in the case when the electron wavelength is much larger than the 
radius of the potential the scattering cross-section is proportional to the electron wavelength. This 
leads to the mobility independent on the electron concentration. In contrast with the case of single- 
layer, neutral and charged defects are, in general, equally relevant for the resistivity of the bilayer 
graphene. 
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Bilayer graphene, that is a two-dimensional allotrope 
of carbon formed by two graphite atomic sheets lli, is a 
subject of hot interest now 0,11, i, d, i, 0, i, i, E,0 El 
motivated by anomalous character of the quantum Hall 
effect 0, B and electron transmission through potential 
barriers [3| due to electron chirality and the Berry phase 
27r, possible use of the bilayer graphene as a tunable-gap 
semiconductor •6] and its other unusual physical proper- 
ties (for review, see RefjlCl). At the same time, it is less 
studied than the single-layer graphene (isl . IT3 | . In partic- 
ular, almost nothing is known about mechanisms of scat- 
tering determining the electron transport in the bilayer 
graphene. Here we consider this problem theoretically. It 
will be shown that for any kind of point defects with small 
enough concentration their contribution in the resistivity 
is inversely proportional to the charge carrier concentra- 
tion resulting in the concentration-independent electron 
mobility. In a framework of perturbation theory, this re- 
sult has been obtained earlier in Refi( see their Eq.(52)); 
we have generalized it on a case of a strong impurity po- 
tential. This situation is essentially different from the 
single-layer case when the scattering by Coulomb po- 
tential of charge impurities leads to the concentration- 
independent mobility whereas the short-range scattering 
centers are almost irrelevant H, H, 17, 18| . 

The bilayer graphene in a simplest approximation can 
be considered as a zero-gap semiconductor with parabolic 
touching of the electron and hole bands described by the 



single-particle Hamiltonian 
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where pi = —ihd/dxi are electron momenta operators 
and m ~ 0.054me is the effective mass, rrie being the free- 
electron mass. This description is accurate at the energy 
scale larger than few meV, otherwise a more complicated 
picture including trigonal warping takes place; we will re- 
strict ourselves only by the case of not too small doping 
when the approximate Hamiltonian ([T]) works. Two com- 
ponents of the wave function are originated from crys- 
tallographic structure of graphite sheets with two car- 
bon atoms in the sheet per elementary cell. There are 



two touching points per Brillouin zone, K and K' . For 
smooth enough external potential, no Umklapp processes 
between these points are allowed and thus they can be 
considered independently. 

The Fourier component of the impurity potential with 
dimensionless charge Z at small enough wave vector 
equals 



V{q) = 



2TTZe^ 
e(q + k)' 



(2) 



where k = 2'iTe'^N (Ep) /e is the inverse screening ra- 
dius, e ~ 2.5 is the dielectric constant due to quartz sub- 
strate, and N (Ep) is the density of states at the Fermi 



energy Ep 



In the model ([T]) 



where we take into account contributions from two spin 
projections and two valleys. Due to the smallness of the 
effective mass the screening radius is 4.5 times larger than 
the nearest-neighbor interatomic distances which makes 
the single valley approximation accurate enough. At the 
same time, for any reasonable doping the Fermi wave vec- 
tor kp <^ K so one can assume that the electron wave- 
length is much larger than the scattering potential radius. 

Let us consider the case of small concentration of point 
defects (to be specific, we will call them impurities) with 
the concentration Uimp and the angle-dependent scatter- 
ing cross-section a {(/)) . Then the defect contribution to 
the resistivity p reads [l^, [20| 



e^vlN{Ep)T{kF)' 



(kp) 



, da (0) 



(1-COS0) (3) 



where vp — hkp/m is the Fermi velocity, t is the mean- 
free-path time. Note that the product vpN (Ep) is pro- 
portional to kp — \/Tm (n is the electron concentration) 
for both single-layer and bilayer graphene, as well as 
for conventional two-dimensional electron gas and thus 
any essential difference in their transport properties can 
be related only to the behavior of the scattering cross- 
section. 
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The expression ([3]) is derived from the standard Boltz- 
mann equation and does not take into account locahza- 
tion (or antilocahzation) corrections which can change 
the results drastically in the regime of small doping when 
the resistivity is of order oih/e^ 21, 2^, 23, 24| (recently 
the problem of weak localization has been considered also 
for the bilayer graphene We will restrict ourselves 

only by the case p ^ h/e^; formal derivation of the Boltz- 
mann equation f or g raphene in this regime will be pub- 
lished elsewhere [25|- 

To determine the scattering cross section one has to 
solve the two-dimensional Schrodinger equation with the 
Hamiltonian ([T]) plus impurity potential V (r) which, af- 
ter simple manipulations (cf. RefQ) can be written in 
the form 



d_ 

dr r 
d l + l 
dr r 



d I 
dr r 
d 1 + 2 
dr r 



91 



ft2 

2mV 



91- 



(4) 



where I = 0,±1,... is the angular-momentum quantum 
number, gi (r) e*'"^ and /; (r) e'('+2)"* are components of 
the pseudospinor wave function, r and (f> are polar coordi- 
nates; to be specific we will consider the case of electrons 
E = h^k^/2m > 0. 

Modifying a standard scattering theory [2^ for the 
two-dimensional case one should try the solutions of 
Eq.Q outside the region of action of the potential in 
the form 



91 (r) 
fiir) 



A 
A 



Ji{kr)+tiHl^Ukr) + ciKi{kr) 
Ji+2 (kr) + UhIH (kr) + ciKi+2 (kr) 



(5) 



where the terms proportional to Bessel (Hankel) func- 
tions describe incident (scattering) waves; the terms pro- 
portional to the Macdonald functions are analogous to 
the exponentially decaying solutions in the case of po- 
tential barrier [3|]. To calculate the scattering cross sec- 
tion one has to find the current operator j and 
its normal component j„ = cos + jy sin (j). The result 
reads: 



Jn = 



hk 

m 







-2i</) 





(6) 



The Bessel and Hankel functions in Eq.([5]) correspond 
to the expansion of the incident plane wave and scattered 
radial wave, respectively. Calculating the average value 
of the current operator Q over the scattered wave we 
find for the cross section 



da{(l)) 



2 

nk 



^ — — oo 



(7) 



which is formally the same expression as for the case of 
single-layer graphene 



18| 



The Schrodinger equation ([4]) has as important sym- 
metry with respect to replacement / < — > g, I < — > —1 — 2 
which means ti — t_;_2- This is the consequence of chiral 
properties of electrons with the Berry phase 27r; a similar 
identityfor the single-layer case with the Berry phase tt 
reads [iSj t; = t-i-i- Thus, Eq.([7|) can be rewritten in 
the form 



da{(t)) 
d4> 



2 

Trk 



oc 



t_i+2> COS [(/ -t- 1) , 



(=0 



(8) 



To understand the behavior of the scattering parame- 
ters ti (fc) in the interesting limit fc — > one can consider 
the simplest case of the potential V (r) = Vq at r < a 
and V (r) = at r > a. Strictly speaking, a sharp jump 
of the potential with atomic scale is beyond applicability 
of our approach since it will induce Umklapp processes 
between the valleys. We assume that the boundary is 
smooth enough in comparison with the interatomic dis- 
tance but much thinner than the electron wavelength (cf. 
Ref.f3). The solution outside the potential well has the 
form ([5]), with A = I and the solution for r < a regular 
at r = can be tried as 



gi (r) = aiJi (qr) + PJi (qr) , 

/; (r) = cr [aiJi+2 (qr) + l3iIi+2 {qr)] 



(9) 



where a = sign {E — Vq) and q — ^2m |i? — Vo|/ft is the 
wave vector inside the well. Using boundary conditions of 
continuity of the wave functions and their first derivatives 
at r = a one can find the scattering parameters ti as well 
as c; , ai and (cf . the case of one-dimensional potential 
i). 

For the case I = —1 taking into account iden- 
tities Ki{z) = K^i{z),Ii{z) = /_i(z),Ji(z) = 



r(i) 



= — (z) one can prove im- 

\2 



— J-i (z) ,and Hi 

mediately that c_i — and i_i cx (fca)^ at /ca — *■ so 
as we will see this contribution in the scattering cross 
section is negligible. Using asymptotic of the Macdonald 
and Hankel functions for I > 2, z —^ 



Ki{z) 



Hi'' (-) 



1/2 
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2)!, 

(1-2)1, 

(10) 



one can prove that for I > 1 and fca — > both ti and c; 
are, at least, of order of (fca)^'or smaller and thus only s- 
channel {I — 0) contributes in the scattering cross section 
so that Eq.® can be rewritten as 



da (0) 
~df~ 



Trk 



|io(fc)|^cos^ 



(11) 
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For the single-layer graphene, <j {(j>) cx cos^ cj)/ 2 and the 
back-scattering is forbidden. On the contrary, for the 
case of bilayer there is a strong suppression of the scat- 
tering at ~ 7t/2. 

For the case I — the wavefunctions outside the well 
([5]) has the asymptotic form 



gi (r) = l + to+ To 



fdr) 



hiy +7) + O {{krflnkr 




O (^{krflnkr 



(12) 



where 7 ~ 0.577... is the Euler constant, tq — — cq. 
Using this we find that to (fc) tends to a finite com- 
plex number ^jig (fc)|^ < 1^ at A; ^ 0. Substituting 

this into Eqs. lfTTjl and ([3]) one can find an estimation 
for the resistivity p ~ [h/ie^) mmp/n. It seems to be 
in a qualitative agreement with the dependence of the 
resistance of bilayer graphene on the gate voltage mea- 
sured in Rcf, llj. The same dependence of the resis- 
tivity on the charge carrier concentration takes place 
for the single-layer graphene with Coulomb scattering 
centers whereas the point defects with short-range 
potential give much smaller resistivity 18| of order of 
p ~ (/i/4e^) nimpO^ ■ For the case of bilayer, on the con- 
trary, there is, in general, no essential difference between 
charge impurities and neutral point defects such as, say, 
atomic-scale roughness of the substrate. 

It is interesting to mention that the scattering by the 
short-range potential in the case of bilayer graphene is 
more efficient than not only in the case of the single- 
layer graphene but also for the conventional nonrelativis- 
tic two-dimensional electron gas where to (fc) oc 1/ In (fca) 
at fca ^ and thus [H [13, [ill 



4e2 nln^ {kFo) 



(13) 



To summarize, we have proven that the scattering by 
point defects in bilayer graphene is more efficient than 
both in single-layer graphene and in conventional elec- 
tron gas. The difference with the single-layer case is just 
due to vanishing density of states for the massless Dirac 
fermions whereas for the bilayer graphene (as well as for 
the conventional electron gas) it is constant. However, 
for the two-dimensional nonrelativistic electrons an arbi- 
trary weak potential leads to formation of a bound state 
in the gap [2^ which results in the logarithmic singu- 
larity of the scattering amplitude at small energies (see 
Eq. (fT3)l V In the case of the bilayer, there is no gap and 
thus no localized states. As a result, the resistivity should 
be just inversely proportional to the Fermi energy, or, 
equivalently, to the charge carrier concentrations. This 
seems to be in agreement with the recent experimental 
data EH. 
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